Effect of fluctuations on vortex lattice structural transitions in superconductors 
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The rhombic-to-square transition field Hn{T) for cubic and tetragonal materials in fields along 
[001] is evaluated using the nonlocal London theory with account of thermal vortex fiuctuations. 
Unlike extended Ginzburg-Landau models, our approach shows that the line Hn (T) and the upper 
critical field IIc2{T) do not cross due to strong fluctuations near Hc2{T) which suppress the square 
anisotropy induced by the nonlocality. In increasing fields, this causes re-entrance of the rhombic 
structure in agreement with recent neutron scattering data on borocarbides. 



PACS numbers: 74.60.-w, 74.60.Ge, 74.70.Dd 

Complex behavior of vortex lattices (VLs) even in cu- 
bic superconductors like Nb has been known for a long 
time Q. Recent progress in understanding the evolu- 
tion of VLs with the magnetic field H and temperature 
T became possible due to availability of large high qual- 
ity crystals of borocarbides which triggered a score 
of small-angle neutron scattering (SANS) 1^-0], scanning 
tunneling ||^,|^ , and decoration experiments Of a par- 
ticular interest is the ubiquitous structural transition be- 
tween rhombic and square VLs in increasing fields along 
the c axis of tetragonal borocarbides [^|-^ • 

For this case, the vortex repulsion is isotropic within 
the standard London and Ginzburg-Landau (GL) theo- 
ries, so that vortices should always form the hexagonal 
Abrikosov lattice, which provides the maximum vortex 
spacing for a given flux density B/(l)o {(f)Q is the flux 
quantum). There is no coupling between the VL and the 
crystal in these models; as a consequence, the VL orien- 
tation is arbitrary and no VL structural transitions are 
expected. A full microscopic theory of the mixed state 
contains this coupling, but involves self-consistent calcu- 
lations of the gap and current distributions, a formidable 
task even for materials with the GL parameter k ^ I . 
The situation simplifies in high-K materials (like boro- 
carbides), for which one can utilize a more transparent 
nonlocal London model |10|. Within this approach, the 
VL coupling to the crystal is provided by the basic non- 
local relation between the current density and the vec- 
tor potential, Ja{r) = f Qaf^i^ — r')A^(r')d'^r', where 
the kernel Q depends on the Fermi surface [ p^pl| , the 
pairing symmetry [|l2|, and the field orientation. Here, 
we consider cubic or tetragonal s-wave materials in fields 
along the c axis so that Q{v) has the square symmetry. 

The kernel Q{r) decays over the nonlocality range 
P = /(2^j^)Coi where £ is the mean-free path for non- 
magnetic scatterers and is the BCS zero-T coherence 
length. The function / decreases slowly with T and is 
suppressed strongly by scattering The nonlocal- 



ity adds to the intervortex interaction a short-range po- 
tential V{x, y) with the symmetry of the crystal. In the 
low field limit, V is irrelevant and the VL is triangular; 
still, V removes the orientational degeneracy and locks 
the VL onto certain crystal direction. With decreasing 
intervortex spacing a{B), the potential V drives the tri- 
angular VL into a square at a field Hn (T) . The transition 
curve Hn (T) is the subject of this work. 

The nonlocal London model describes correctly the ob- 
served structure and orientation of rhombic VLs in small 
fields m and the structural evolution toward the square 
Moreover, the fast increase of Ho with increasing 
impurity concentration predicted by the model has been 
verified by SANS [^. Still, there is an open question 
on what happens when the applied field H > Ha keeps 
increasing. The nonlocal London model per se implies 
that the square VL should persist all the way up to Hc2 
since the shorter the intervortex distance, the stronger 
the role of the square-symmetric potential V. This con- 
jecture is supported by calculations based on the GL 
theory extended to include the 4th order gradient terms 
P, p^p^ ; the model predicts that the lines Ho{T) and 
Hc2{T) should cross as sketched in the inset to Fig. 1. 

However, recent SANS data for LuNi2B2C Q indicate 
that instead of crossing Hc2{T), the hue Ho{T) curves 
up to avoid Hc2 and becomes two- valued, see Fig. 1. The 
upper branch of Ha (T) marks the re-entrant transition 
of the rhombic VL in increasing fields, which then evolves 
toward the triangular VL as H ^ Hc2- In this Letter we 
argue that vortex fluctuations combined with nonlocal 
effects can account for this unexpected behavior. 

Below we evaluate the mean-squared amplitude of 
thermal vortex fluctuations using the elastic energy of 
the deformed VL. Unlike the case of isotropic supercon- 
ductors, an important role in this energy is played by 
VL rotations relative to the crystal. We show that 
remains finite at the transition line Hn (T) , but diverges 
as B approaches Hc2{T). In the vicinity of Hc2{T), the 
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anisotropic nonlocal potential V{x, y) is averaged out by 
fluctuations, and the interaction becomes isotropic. As 
a result, the rhombic VL becomes preferable, turning 
into the hexagonal Abrikosov VL as i? — > Hc2- The 
re-entrance of the rhombic VL in increasing fields occurs 
if ~ i.e., the amplitude of fluctuations needed to 
wash out the nonlocal effects is much smaller than that 
required for the VL melting. This re-entrance, therefore, 
can happen in nearly cubic materials, in which vortices 
are not split into "pancakes" and fluctuations are weak. 
We evaluate the shape of the Hu{T) curve using model 
parameters of LuNi2B2C for which in the clean limit 
p{T) ~ (1 ^ 2)^0 0,0. We reproduce the shape of this 
curve seen in the SANS data [p^ . 

Let us consider the equilibrium square VL with diag- 
onals along [100] and [010] directions, which we take as 
X and y axes. We are interested mainly in the high fleld 
region, where the VL can be considered as incompress- 
ible [0. Then the displacement of vortices Ui{x,y\z) 
{i = x,y) satisfy divu = dxUx + dyUy — and the sym- 



metric strain tensor Ui 



{djUi -t- diUj)/2 has only two 



independent components: u^x and Uxy Since the VL ori- 
entation is locked on the crystal, rotations of VL about 
[GDI], LOxy = (dxUy — dyUx)/2, cause an energy increase. 
Then the energy density U of the deformed VL reads [19| : 



2U ^Csul^ + Cxuly + c^ujly 



044(^2 u)^ 



(1) 



Our notation is motivated by the following. A uniform 
displacement = ex, Uy — —ey (e ^ 1 is a constant) 
is of a special interest: it is this deformation that trans- 
forms the square above Hn into a rhombus for B < Ha ■ 
This deformation was named "squash" ||20[] , so that is 
the squash modulus. Since this is a second order phase 
transition, Cg must vanish at B — Hn- Further, the "sim- 
ple" a;-directed shear is Ux = ey, Uy — so that only 
Uxy 7^ 0; thus the notation Cx for this shear mode. The 
shear energy depends on the shear direction, but the cor- 
responding moduli can all be expressed in terms of Cx 
and Cui [|9| . We take all moduli in Eq. (|]) except the tilt 
C44 as practically nondispersive . 

Next, we calculate the mean-squared vortex displace- 



ment u^iT, B) 



Writing the total elastic energy 



in the Fourier space and utilizing the equipartition theo- 
rem, we obtain |pi|: 



f qdqdkz 



2T 



(27r)3 q^c{<t>)/A + CiA{k.,q)kl 



(2) 



where kx — qcoscf), ky ~ qsmcf), and c(0) — Cuj + 
Cx sin^ 2(j) + Cs cos^ 2(j). With the large flelds tilt modu- 
lus C44(fc) = C44(0)/(l -I- A^fc^), integration over kz and q 
gives — UqIJ. Here Uq ~ TX/a^ ^y CxC44{0) is the mean- 
squared displacement for the triangular VL with cgg = Cx 
and C44(0) = B^/Stt ||2|,|i^. Taking for a hexagonal 
VL in a uniaxial material for Bjjc 18|, we arrive at 



m2 ^ 16V2n^X,X,^,TT]/(j)ls{b) , 
=b{l-bflii{2 + l/V2b). 

Here, A^.c sue penetration depths, b — B/Hc2, and 



V 



2^x r 



(3) 
(4) 



(5) 



where K(m) is a complete elliptic integral. 

As follows from Eqs. @)-@,y^{T,B) remains finite 
at the instability point Cs — 0. Thus, vortex fluctuations 
do not affect the mean-field character of the second order 
transition at Hn ■ This unusual behavior occurs because 
the rotational modulus c^^ is finite, otherwise v?' would 
diverge logarithmically as S ^ Hu ■ 

The moduli Cx and in intermediate fields, Hd 
B ^ Hc2, are estimated as Cx — 4>oB / [SttX)"^ and 

. For c^/cx ^ p^/a^ < 1, 



p2s2/27rA2\/6 



we obtain 77 ~ (2/7r) ln(16cx/ctij). Hence, the nonlocality 
affects via the weakly-varying factor 77 ~ ln(^/p6^/^), 
i.e., u'^{T,B) is insensitive to detailed behavior of 
c^{T,B) and Cx{T,B). Since ^ cx (i?c2 - B)-^/"^ di- 
verges only at -ffc2, the fluctuations suppress the weak 
VL-crystal coupling mainly near Hc2- 

We now turn to evaluation of the curve Ha{T). The 
free energy F of fluctuating VL can be written as: 



F ■ 



B^ 
Stt 



-(CC^+«2)gV2 



X-^G^ + P^X^G^Gl 



(6) 



Here the sum runs over Gx — 7r\/2(m — n)/atan-'^/2(/3/2), 
Gy — TT^/2{m + n) tan^^^ {(3 /2) / a, which form the recip- 
rocal lattice; /3 is the apex angle of the rhombic cell, 
and TO, n are integers. The last term in the denominator 
describes the nonlocal correction to the London theory 
0,0. The factor exp(-C£2GV2) accounts for the fi- 
nite size of the vortex core |18| ] , which provides a cutoff in 
the London model. The B dependent quantity C was es- 
timated theoretically as 0.5 < C < 4 |^J]. By and large, 
the data on the field dependence of SANS form-factors 
are consistent with C ~ 1 [0| , the value we adopt below. 

Fluctuations enter the energy (^) via the Debye- Waller 
factor exp{—u'^G^/2), which accounts for the thermal 
smearing of the vibrating VL |2^. For > p^, the 
fluctuations wash out the nonlocal corrections in F, mak- 
ing the averaged vortex interaction isotropic. Hence, the 
re-entrant transition can be driven by thermal vortex dis- 
placements of the order of p ~ ^o- The instability fleld 
Ha (T) thus can occur well below the melting line, which 
in materials of interest here is rather close to Hc2- 

To find Ha{T) at which = d^F/dp^\p^^/2 = 0, we 
differentiate F of Eq. (p|) and obtain: 



^^[(2pTOn)^ + (^ 



0. (7) 



2 



Here, p = [nC + xiT)r){t,b)/ s{b)]b, g = m'^ + n^, d = 
fj, + g + Cib){m'^ — n^)^, fi = l/27r6K^. The dimension- 
less control parameters x ^-nd C quantify the amplitude 
of thermal displacements and the nonlocal corrections: 



X 



2 U 



(8) 



Note that nonlocality enters also the exponent p in Eq. 
via the parameter 77. 
For the further analysis, we assume that A(T) = 
A(0)/(1 - t2)i/2 and ^{T) = ^(0)7(1 - t^^^, where 
t = T/Tc (qualitatively, our results do not change if other 
plausible T dependences are used). Then 



16^/27r3A,(0)A6(0)T, 



Xo 



C = Cob{l-t^), (9) 



In the clean limit Co ^ 1 ; with increasing scattering, Co (i) 
drops fast 0. For LuNisBsC with = 16 K, « 70A, 
Aa(0) « lO^A, and Ac(0) w 1.2 x lO^A, we estimate 
Xo ~ 6.4 X 10~^ <C 1. The smallness of xo indicates 
that fluctuations contribute little to the thermodynam- 
ics of stable VLs ||2^,|l^. As shown below, being crucial 
on the upper branch of Ha (T) , fluctuations are negligible 
on the lower branch. 

Strictly speaking, Hn{T) should be calculated self- 
consistently taking into account the effect of fluctuations 
on the relevant moduli [ p6[ . However, since is finite 
at Hn, we may neglect the thermal softening of Cx and 
c^. Then, Hn(T) is just a root of Eqs. (|^), which we 
find numerically. The factor 77 which enters p in Eq. (|^) 
is a much weaker function of t and b than x{t)/s{b) {rj 
varies from 1.6 for c^j/cx = 1 to 3 for c^j/cx = 0.1). For 
this reason we disregard variation of 77, adopting rj = 2.7 
for LuNiaBzC with c^/cx =i 0.2 a.t Hn 

The results of the numerical solution of Eq. (0) are 
shown in Fig. 2. It is seen that fluctuations do give rise 
to the re-entrant square-to-rhombus transition in high 
fields, in a qualitative agreement with SANS data of Fig. 
1. In fact, fluctuations change radically the VL phase 
diagram in high fields, while weakly affecting the low- 
field branch of Ha- The difference between Hn{T) and 
Hc2{T) is significant (except the low-T clean limit) which 
justifies the use of the London model. As the ratio po/Co 
decreases (e.g., due to nonmagnetic impurities), the re- 
gion of the square VL on the H — T diagram shrinks. 
The raise of the lower branch of Hn (T) has been seen on 
Lu(CoxNii_x)2B2C, for which the mean-free path £ was 
suppressed by the Co doping ||l3| . 

It is worth noting that although the calculated curves 
Hn {T, t) reproduce correctly qualitative features of the 
SANS data of Fig. 1, actual position of the upper branch 
is sensitive not only to the accuracy with which we know 
the elastic moduli and the parameters for their evalua- 
tion, but also to the precise value of the empirical cutoff 



constant C, see Eq. (|i|). The information on the actual 
position of the upper branch of Hn {T, £) is still scarce, 
and we hope to refine our approximations when the data 
are available. 

Now we comment bricfiy on the possible effect of the 
VL transition on the flux pinning. Since the instability of 
the square VL at Hn is not accompanied by divergence 
of u^, the critical current density Jc evaluated within the 
collective pinning theory, should not be sensitive to the 
VL transition. Indeed, the correlation function of vortex 
displacements {u{r)u{r')) can be evaluated with the help 
of Eq. (H) in which T is replaced by_7p exp 7k(r — r'), 
where 7^ is the pinning parameter p2| . At Hn, the 
squash modulus vanishes, but {u(r)u{r')) remains of the 
same order as for a triangular London VL, to the accu- 
racy of the weak logarithmic factor 77 ~ ln(^/p6^/^) ~ 1. 
Thus, contrary to the claim of Ref. |l^, neither the pin- 
ning correlation length nor Jc are significantly affected 
by the squash softening near Hn ■ 

The two- valued Hn{T) and the re-appearance of the 
triangular VL at — > Hc2 are generic features which 
are not limited to nonmagnetic borocarbides. The low- 
T SANS experiments on antiferromagnetic TmNi2B2C 
have revealed the triangular VL near Hc2 , which evolves 
into a square as the field decreases pTf . At this stage 
the effect of antiferromagnetic ordering upon VLs is un- 
clear, but ffuctuations could certainly contribute to the 
re-entrant VL transition in TmNi2B2C as they do in 
LuNi2B2C. Similar behavior was seen in YNi2B2C [ p8[ . 
Another candidate for studying effects of vortex fluctu- 
ation is VsSi, in which the rhombic VL was observed at 
T < 5°K and H = lOkOe {Hn ^ 15kOe); as T increases 
at the fixed field, the rhombic VL evolves toward the 
hexagonal one as T ^ Tc2{H) [Q. 

In conclusion, we present a model of the structural 
VL transition at Hn{T) affected by thermal fluctua- 
tions of vortices. We show that the curves Hn{T) and 
Hc2{T) do not cross, instead Hn{T) becomes two- valued. 
The upper branch of Hn [T) corresponds to a re-entrant 
transition of the rhombic VL, in accordance with recent 
SANS observations on LuNi2B2C. 
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FIG. 1. The transition line Hn(T) (circles) in LuNi2B2C 
observed by SANS [17]. The inset shows Hn{T) (dashed line) 
predicted by the extended GL theory [6] without vortex fiuc- 
tuations. The solid lines shows Hc2{T). 




FIG. 2. The transition lines Ha (T) obtained by numeri- 
cally solving Eq. (^) for xo = 0.0064, C = 1, and a few ratios 
of p/^o- The dashed line is Hc2(J). 
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